Analytic Aperture Calculation and Scaling Laws for Radio Detection of 

Lunar-Target UHE Neutrinos 
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ABSTRACT 

We derive analytic expressions, and approximate them in closed form, for the effec- 
tive detection aperture for Cerenkov radio emission from ultra-high-energy neutrinos 
striking the Moon. The resulting apertures are in good agreement with recent Monte 



Carlo simulations and support the conclusion of James Sz Protheroe (20091 that neu- 



trino flux upper limits derived from the GLUE search (Gorham et al. 2004 1 were too 



low by an order of magnitude. We also use our analytic expressions to derive scaling 
laws for the aperture as a function of observational and lunar parameters. We find 
that at low frequencies downward-directed neutrinos always dominate, but at higher 
frequencies, the contribution from upward-directed neutrinos becomes increasingly im- 
portant, especially at lower neutrino energies. Detecting neutrinos from Earth near the 
GZK regime will likely require radio telescope arrays with extremely large collecting 
area {Ae ~ 10^ m^) and hundreds of hours exposure time. Higher energy neutrinos 
are most easily detected using lower frequencies. Lunar surface roughness is a decisive 
factor for obtaining detections at higher frequencies {i^ ^ 300 MHz) and higher energies 
(£;> 10^1 eV). 



1. Introduction 

The ubiquitous presence of ultra-high-energy (UHE, > 10^^ eV) cosmic rays suggests the 
existence of an equally ubiquitous and similarly high-energy cosmic neutrino population, either as 
a result of the various mechanisms for generating muon neutrinos from charged pion decay in the 
vicinity of the cosmic-ray acceleration region ( Bahcall Sz Waxman|2001 ), or during interactions with 
the cosmic background radiation (GZK effect, Greisen||1966 Zatsepin Sz Kuz'min l|l966l ). One effort 
to detect such neutrinos involves radio observations of the expected Cerenkov burst emission when 
these neutrinos interact with the Moon. Several experiments have already attempted to detect this 



signal (e.g., Hankins et al. 1996; Gorham et al. 2004; Beresnyak et al. 2005 Buitink et al. 20081 
all with null results to date. But even null results can be converted into useful upper bounds for 
the UHE neutrino flux, provided that the effective aperture, the area times the solid angle through 
which incident neutrinos are detectable, is known. 



For future experiments, it is clear that the larger the aperture, the greater will be the possibility 
of achieving a detection, or the more decisively constraining will be the inferred upper limit. Of 
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particular note is the fact that all the experiments so far have suffered from particularly weak 
coverage of the energy domain near the GZK cutoff near lO^^'® eV, the region of perhaps the 
highest cosmological interest, and also possibly a local peak in the neutrino spectrum. The poor 
coverage is largely due to inherent limitations on radio observations, but it may also in part be due 
to the difficulty in achieving optimal tailoring of the observing parameters. Having access to closed- 
form expressions of general validity would yield not only a convenient means for calculating the 
aperture, it would also assist in such optimization of experimental design, for penetrating deeper 
into this hitherto uncharted neutrino regime. 

Despite their potential value, no such closed- form expressions are currently available in the 
literature, and so the primary purpose of this paper is to provide such expressions. We then 
apply them to the question of how to optimize for detection of GZK and other UHE neutrinos 
by looking for radio Cerenkov signals from the Moon. The approach is general enough for future 
modifications to accommodate other types of experiments, such as terrestrial ice sheets seen from 
airborne balloons, or the Moon seen from lunar orbit. 

To account for the neutrino properties and their detectability in closed-form aperture expres- 
sions, we rely on previous determinations of the basic attributes of radio Cerenkov emission from 
ultra- high-energy (UHE) hadronic showers (e.g., James &: Protheroe|2009 ) and reviews of the basic 
physics of the Askar'yan effect, whereby neutrinos generate charge excesses in hadronic showers 
(e.g., Alvarez-Muniz &; Zas 2001 1 . Here we will not comment on these physical processes, but 



merely quote the results from the literature as we apply them to our specific problem. 

In this paper we assume that all contributing neutrino showers occur in the lunar regolith 
near the surface, with a fixed index of refraction of rir = 1.73 (James &: Protheroe 20091, so we 



ignore inhomogeneities such as the sub-regolith. James Sz Protheroe (20091 have evaluated a sub- 



regolith contribution and find that it can only significantly increase the aperture for Earthbound 
detection at very high neutrino energies {E > 10^^ eV), because only such high-energy neutrinos 
produce sufficiently strong fields to be above the minimum detectable level at the surface. Of 



course, experiments from a closer distance, such as lunar orbiters like LORD (Gusev et al. 20061 



could receive significant contributions from the sub-regolith, and ignoring gradients in the lunar 
material introduces potential errors but is certainly a greatly simplifying assumption. 



2. General form of the aperture calculation 

Our goal is to determine the rate of detection of hadronic showers initiated in the Moon by the 
capture of a high-energy cosmic neutrino, via the detection of the resulting radio-frequency electric 
fields propagating from the Moon. One way to conceptualize this detection rate is in terms of an 
aperture size, which when multiplied by the incident neutrino flux in each energy bin of interest, 
gives the detection rate in that bin. Since the diffuse neutrino background flux is scaled per area 
and per solid angle, the aperture will be in units of an area times a solid angle (Williams 2004 1 . 
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Our approach for analytically specifying this aperture begins with identifying the maximum 

possible aperture A^ax that an amount of mass equivalent to the lunar mass could possibly achieve, 
assuming that each neutrino interacting with that material is detected exactly once. Note this 
would actually require the lunar material be unphysically spread out, among other impossible 
requirements. Thus it is merely a starting point, and it follows clearly that 

Amax{E) = At,Mk{E) (1) 

for mass M (of the Moon) exhibiting a cross section per gram k{E) for initiating charged particle 
showers for neutrinos at energy E, incident over the full 47r steradians of sky. 

Our approach is then to reduce the aperture by eliminating events that are blocked from 
occuring by virtue of prior absorption of the neutrino elsewhere in the spherical Moon, and then 
reduce it further by requiring that the events yield radio signals above the detection threshold 
of the specified telescope system along some ray that intersects the detector. These reductions 
are substantial for three reasons: ultra-high-energy (UHE) neutrinos are drastically truncated by 
the opaqueness of the Moon, total internal reflection at the lunar surface reduces the rays that 
successfully cross the falling refractive index at this boundary, and rays from depth in the Moon 
are signiflcantly attenuated by lunar radio absorption. The remainder of this paper is devoted to 
quantifying the aperture reductions stemming from these three effects. 



2.1. The phase-space psirtition 

In order to include these corrections, we first subdivide Amax by multiplying it by the frac- 
tional (normalized to unity) phase-space volumes of all the processes that contribute to Amax, and 
then weight each subdivision by an efficiency fraction, or probability, by which that phase-space 
component contributes to the observable aperture. Hence, the detection rate for a neutrino flux 
distribution I{E, u), where I{E, u) is per energy bin at energy E and per solid angle along direction 
v and per target area, is 

n = JdVmJ dCl^ j dh^J dE A^ax{E) I{E,U)e-^^ Hd Hr , (2) 

where the phase space consists of three dimensions of lunar volume dVm that account for all possible 
shower locations inside the Moon, two dimensions in dO-y that account for the possible directions 
of the rays along which the electric fleld of the shower can propagate (after leaving the Moon), 
and two dimensions in dCt^,, which account for the possible directions of the incident neutrino. 
The weighting function e""^" accounts for the penetrating fraction of neutrinos which reach the 
phase-space contribution element under consideration (the first correction mentioned above), the 
Heaviside step function Hr selects the outgoing radial rays from d(l^ that do not totally internally 
reflect at the lunar surface (the second correction), and Hd further selects the rays from dA-y that 
are bright enough to detect (the third correction) . 
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2.2. Surface roughness effects 

To determine whether or not a particular ray will escape total internal reflection at the lunar 
surface, it is necessary to know the angle that the ray meets the surface, which can be altered by 



lunar surface roughness on the wavelength scale or larger. Shepard et al. (1995) analyzed radar 
reflections from the lunar surface and found that the surface irregularities are self-similar (fractal), 
and the root-mean-square roughness angle, which we convert to a Gaussian halfwidth by multiplying 
by \/2, can be parameterized as 

ao{X) = V2 tan"^ (0.29 A'^'^^) = ^/2 tan-^(0.14 u^-^^) (3) 

where A is the spatial scale in cm, in GHz, and cjo is in radians. 

Even though the surface may be tilted symmetrically in any direction, the way surface rough- 
ness alters the local refraction angle may have an important impact because of the nonlinearity 
of the Heaviside step functions. This nonlinearity can be particularly important for "downward" 
neutrinos, which are the most likely to produce total internally reflected rays becasue of their 
downward-pointing Cerenkov cones. To illustrate this. Fig. [T] depicts schematic ray paths of the 
escaping Cerenkov cone both with and without surface roughness. In general, favorable surface 
tilts increase the number of rays escaping from downward-directed neutrinos more than the unfa- 
vorable tilts reduce it, because there may already have been little or no signal prior to the inclusion 
of surface tilt (and note that once a signal contribution becomes zero it never becomes negative, 
so unfavorable tilts come with little penalty). The enhancement becomes especially important as 
the root-mean-square surface roughness angle exceeds the Cerenkov width, i.e., at high observing 



frequencies (z^ ^ 300 MHz, see section 4.1 ). On the other hand, we show below that whenever there 
is significant contribution from "upward" neutrinos (neutrinos that have survived a significant se- 
cant of lunar rock and approach the surface from an upward angle) , surface roughness plays a less 
important role. 

Accounting for surface roughness may be accomplished by using a probabilistically smoothed 
TCr function defined by 

/tt/2 /-co /"''"/2 foo 

dcp' / dw e-'^'n^,^^-^^ = 27r-'/^ / / dw e-'^'w^,^, , (4) 

which allows the surface normal to be tilted, relative to the radio ray incident along 7, by a small 
polar angle a = ctoW in a random azimuthal (</>') direction. The polar angle a is approximately 
normally distributed with Gaussian halfwidth Uo, as supported by the findings of [Shepard et al.] 



(19951 in the limit of small a. The expression (fio- • 7)/(n • 7) accounts for the higher likelihood 
of encountering a tilt that reduces the angle of incidence than one that increases it, owing to the 
relative increase in projected area of the former, but this efl'ect is deemed to be too negligible to 
track when the tilt angles are small, consistent with our other approximations. Here 'Hw,(j>' expresses 
the no-internal-reflection selection rule in the form of a step function whose argument is positive 
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Fig. 1. — (Left panel) Side and face-on views of Cerenkov cone ray paths escaping the lunar surface 
for downward-directed neutrinos, assuming a smooth surface. (Middle) Same, but including surface 
roughness. Note that a larger fraction of the Cerenkov cone escapes. (Right) Same, but for upward- 
going neutrinos. 



whenever the values of w and (j)' allow the ray to cross the surface, so TCr itself becomes a smoothed 
function rather than a formal Heaviside function. 

It should be noted that we choose the somewhat nonstandard convention of allowing w to 
be either positive or negative, corresponding to surface tilts that either help or hinder ray escape, 
and take the azimuth of the surface normal to wrap only from — 7r/2 to tt/2 (further applying the 
left/right symmetry to restrict to one azimuthal quadrant from to it/2), rather than restricting the 
tilt angle to positive values and wrapping around the full 27r in azimuth. We choose this convention 
because it indicates more explicitly when the tilt is helping or hindering escape, an issue to which 
the Heaviside functions are highly sensitive. Also note that Tin is normalized to unity in the purely 
hypothetical situation where Tl'aj,4>' always unity, as desired. 

Including surface roughness can increase the aperture significantly, but it also adds two addi- 
tional dimensions to the integral over fractional phase space, and it introduces a new angular scale, 
(To, to the several other important angular scales that appear in the calculation. It is natural to 
expect this complication to be worthwhile whenever cto is appreciable relative to the other angular 
widths that affect the aperture, for example when observing at high frequencies where the Cerenkov 
cone is thin. The scaling laws we derive bear out this expectation. 
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2.3. Isotropic neutrino flux 

As the current interest is primarily in cosmologically distributed sources, rather than targeted 
sources, we focus on isotropic neutrino flux distributions I{E), and simphfy the form of the aperture. 
Then we may write the detection rate TZ as 



n = J dE I{E)A{E) , (5) 
where the energy-dependent isotropic aperture is given by 

A{E) = Amax [ dVm I dCl^ I dQ, e-"'^ Hr Hd ■ (6) 



Again note that each fractional phase space integral is by definition normalized to unity, and 
^max specifies the units of A{E), but not the scale of A{E). The great majority of the fractional 
phase space will generally not contribute to neutrino detections, and thus e""^", TLr, and TLd exert 
stringent constraints on A{E) in ways we will now calculate. 



2.4. Spherical symmetry simpliflcations 

The fractional phase space in eq. ^ includes 7 total dimensions, or 9 including a surface- 
normal tilt distribution, so would certainly represent a daunting undertaking to compute in closed 
form. However, we receive welcome assistance from the basic spherical symmetry of the Moon, in 
concert with an isotropic incident neutrino flux. Capitalizing on that symmetry actually allows us 
to eliminate three of the dimensions of the phase space (two because the lunar aperture operates the 
same way as seen by observers from all directions, and one because the lunar aperture seen by those 
observers has an axial symmetry around the Moon center). This reduces the phase space (even 
with surface roughness) to 6 dimensions, which is certainly more tractable, though still requiring 
extensive use of approximations to achieve a closed- form result. 



2.5. Aperture calculation from the lunar volume-centered perspective 

Several choices are possible for coordinatizing the phase space, each with its various compu- 
tational advantages and challenges. We adopt a Moon-centered perspective that scans in spherical 
coordinates over the volume of the Moon by accounting only for the distance h below the lunar 
surface along the radial direction, and we coordinatize the incident neutrino direction in terms of 
the glancing angle a that the neutrinos make to the surface (so a is the complement of the angle 
of incidence to the normal), so chosen because it tends to be a small quantity. Our convention 
is a > for upward neutrinos that must penetrate a significant lunar secant before reaching the 
detectable zone. 
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Hence, the 4-dimensional configuration phase space is divided into one dimension to describe 
the location of interest in the Moon, one dimension to describe the incident neutrino angle, and 
two dimensions to describe the outgoing electric field rays. At this point in the calculation, we 
temporarily consider the two dimensions of outgoing radio rays to be outside the Moon, but shortly 
we will convert to a coordinatization where these are inside the Moon and along the Cerenkov cone 
of the shower. As mentioned above, in addition to these four dimensions, there are generally two 
more integrations to account for the local surface normal variations on the scale of the assumed 
surface roughness. 

The crucial simplification stemming from the spherical symmetry and the isotropic neutrino 
assumption is that fractional phase-space component weights are based only on the probability 
that a hypothetical observer at Earth distance randomly chosen over the solid angle of the lunar 
sky could detect the radio signal, without specifying any particular viewing angle. This succeeds 
because we are in effect fixing the location where the shower occurs along an arbitrary radial ray 
in the Moon, randomizing the orientation of the observer, and asking what is the probability that 
such a randomly located observer could detect that shower. 



Translating the above into an expression for the aperture yields 

1 r/' , /• ,A 3 /-^ 

-tt/2 



1 /"'^/^ /■ - 3 

A{E) = Amax{E) X - da COS a X I dQ^ x dr e-^-'Ho , (7) 

^ J-Tr/2 J Jo 



where Rm is the radius of the Moon. Note how each fractional phase space component is still 
normalized to unity, in keeping with the partition of Amax over its contributing fractional phase 
space. Defining the neutrino interaction length to be Ly{E) = l/pn^E) for cross section per gram 
k{E) and mass density p (and we take p = 1.8 g cm^'^ in the regolith where the showers occur), we 
obtain from eq. ([T]) 

AmaxiE) = AnMKiE) = Att (^RIp) k{E) = \ ^A^ , (8) 



3 ^"V ' ' 3 
where 

Ao = Att'^RJ (9) 

is the maximum attainable aperture for a spherical object for which neutrinos can be detected 
at most once (see the Appendix for the reason that this "maximum" aperture could actually be 
exceeded by added contributions from neutrinos with initial energies above E). Note that Aq is 



the geometric lunar cross section iiRm^ times the full 47r steradians of illumination, as in Williams 



(20041, and eq. ([7| yields A{E) = Aq if we set TId'Hr = 1 and treat in the highly opaque 



limit, where there is as yet no correction for downgrading of higher-energy neutrinos (again, see 
the Appendix for how the effective is reduced by downgrading) . 

Since we treat only neutrino energies for which the Moon is highly opaque, we have Ao <C Amax{E), 
owing to the high shielding of the lunar interior, and thus Aq makes for a better fiducial reference 
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to use in our aperture expressions. Thus we will express the aperture A{E) in the form 

A{E) = Ao P{E) , (10) 

where P{E) is to be interpreted as the fraction of neutrinos entering the Moon at energy E that 
will actually be detected by the instrument under consideration, assuming all such neutrinos will 
create showers, and corrected slightly for the detection of downgraded neutrinos originally at higher 
E (see the Appendix). 



2.6. Converting from exterior to interior ray angles 

Since detectability of the radio rays is a crucial issue, the field strengths generated by the 
hadronic showers appear prominently in the aperture calculation. However, the characteristics of 
these fields trace the Cerenkov cone inside the lunar material, whereas the ray-angle phase space 
is referenced to the emergent solid angle outside the Moon where the detector is located. This is 
inconvenient, so we now convert the ray angles to being inside the Moon where the field properties 
are more easily expressed. To effect that change, we must not only convert dO^ to an angular 
integral over an interior solid angle, we must also account for the solid-angle magnification factor 
that appears whenever emerging rays refract across a dropping index of refraction, which here falls 



from its internal value of = 1.73 (James & Protheroe 20091 to unity. This magnification factor 
is (e.g., |Gorham et al.]|2004[ ) 

m = , (11) 

V 1 - sin^ fj 

where (3 is the angle of incidence from the normal of the radio ray as it encounters the lunar surface 
from the inside. 

We coordinatize this interior-ray solid angle using the polar angle A (measured relative to 
the Cerenkov angle 9c so the A of interest are usually quite small), and the azimuthal angle (j) 
around the Cerenkov cone (with the convention that = corresponds to the direction nearest 
to the surface). Hence we replace the normalized solid angle dCl^ by the transformed solid angle 
dA sin(^c + A)d(j) ^{[3)/2tt^ where we take (f) only from to vr because of the left/right local 
symmetry. The interior phase space element has its normalization altered by the ^{(3) factor, but 
the resulting exterior solid angle will not exceed its full unit value because the contributing interior 
solid angle is actually quite small, owing to truncation by TLr as seen below. 

Indeed we expect all integrals to be truncated by the selection rules imposed by TLr and TLd-, 
so as a minor convenience we will extend all finite integration limits to infinity. This step has no 
physical significance and will not alter the outcome of the calculation, it merely removes unnecessary 
emphasis from the arbitrary limits of the integrals, and returns the emphasis to the selection rules 
themselves. Since we are working in the limit where the radio rays are rapidly attenuated in the 
Moon, we assume all detections occur near the surface, so we also replace the integral over r with 
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an integral over z = h/L^^ where h is the depth below the surface and is the electric field 
dissipation length (so is the photon mean- free-path) . We also approximate by in the 

z phase-space integration. These are all excellent approximations in the domain of interest of our 
calculation. 

Taking the definition of P{E) from eq. (10), we combine eqs. ^ and (m to yield 

1 L /"oo poo poo 

P{E) = - -1 da cos a dAsm{ec + A) dcp dz e~^'' Hr Tin C ■ (12) 

^ J^i^ J-oo J~oo Jo Jo 

Note that already the scale of A{E) has been reduced by the factor L^/Ly <^ 1, the proportion by 
which neutrinos overpenetrate to depths beyond where detectable electric fields can emerge. The 
remaining factors that will further reduce the aperture, mediated by the TC truncations, will be 
included next, and will invoke additional approximations. 



2.7. Imposing the Hr and Tin selection rules 



The constraint described by 'Hw,(j)' is that the radio ray emergent from the electron shower must 
not internally reflect at the surface of the Moon, and the constraint described by is that the 
signal be detectable by the instrument of interest. Each Heaviside function uses whichever quantity 
must be positive to apply the constraint, which appropriately truncates the limits of integration. 

We first consider the requirement that we count the contribution only from rays that are 
associated with detectably strong radio waves. We assume that each shower induces a Cerenkov 
cone with a field strength along each ray that depends on the deviation A in polar angle from the 
Cerenkov peak angle 9c- If the field strength £ is distributed over A in an approximately Gaussian 



way (small deviations from this are discussed by Scholten et al. 2006 Gusev et al. 2006), then 



£oi\\{P) e 



-{A/A„)2 



(13) 



where Aq is the Gaussian halfwidth of the angular distribution around the Cerenkov angle, £o is the 
strength of the field at the shower along the Cerenkov angle, is the number of radio dissipation 
lengths the field passes through before exiting the Moon at the frequency in question, and t||(/3) 
is the field transmission coefficient appropriate for the diverging rays that fill the emergent solid 
angle of interest. We assume the polarization is in the plane of incidence (termed "pokey" electric 
polarization), as that is the dominant polarization for the rays most likely to escape total internal 
reflection, and f3 is the angle of incidence to the surface normal (inside the Moon). 

The result begins with the standard expression for the field transmission coefficient for plane 
waves ( |Williams|[2004l ) 



/ Ur COS (3 

COS f3o 



(14) 
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where ry is the field reflection coefficient for "pokey" electric polarization, and /3o is the angle of 
refraction relative to the normal (outside the Moon) as the rays pass through the surface of the 
Moon into free space, so 

j3o = sin~"^(nr sin/3) . (15) 



In an appendix of Williams ( 2004 ) , a ray-tracing technique is described for converting this plane- 
wave transmission coefficient to one appropriate for the diverging rays that fill the observable 
emergent solid angle, which agrees with an analytic result she quotes from a private communication 
with Dave Seckel. The result is 



/ tan/3 
tan /3o 



1 - 



.2^ _ 



K = 



where 



tan (3 ^ 2 sin /3 cos /3o 

tan/3o " sin(/?+/3o) cos(/3 — /3o) 

2 sin /3o cos (3 



sin(/3o + 13) cos(/3o - 13) 



(16) 



(17) 



It may be noted that the Williams (2004 1 numerical result, the above analytic expression. 



and the analytic expression cited in Gusev et al. (20061, present sequentially more pessimistic 



transmission coefficients, at roughly the 10% level per step in the sequence. This motivates both our 
using the "intermediate" level of optimism, and the application of the rather crude approximation 
outlined below. It appears that the transmission of diverging radio rays from a hadronic shower 
remains a problem that is not completely solved. 

In this paper we use expressions for the maximum electric field and the Cerenkov cone width in 



the regolith as given by James & Protheroe (20091, which tend to give narrower detection windows 



and lower apertures than previous values in the literature. For the maximum field we have 



£o{d, E) = 0.0845 



V 



m • MHz 



d 
m 



-1 



Es 
EeV 



GHz 



1 + 



2.32 GHz 



-1.23 



(18) 



where d is the distance from the shower (and m is the unit in meters), v is the observing frequency, 
and Eg is the shower energy in EeV (10^^ eV). Approximately 20% of the incident neutrino energy 



is deposited in hadronic showers, independent of neutrino fiavor (James & Protheroe 2009), so we 
take Eg = 0.2E. The Cerenkov cone 1/e half- width is given by 



A, 



0.05 



"GHz" 




V 





(19) 



Note that we have multiplied the angular width constant Cu = 2.4° of (James & Protheroe 



20091 by the factor 1/^/ln 2 = 1.2 to make Aq be the 1/e half- width rather than the half- width 
at half-maximum, as this is more attuned to use with the familiar exponential attenuation factors. 
We also convert the angle units to radians for use in the scaling laws that follow, since radians are 
the useful unit for testing the small-angle approximations that will be invoked shortly. 
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For a telescope with a field detection threshold Smin, the selection rule to have a detectable 



signal in eq. (13) is 

Srmn < f o^n (/5)e-(^/^°)' e"^- , (20) 



or, solving for A, 

A < Ao 



where s is the path length from the shower to the surface and is the electric field dissipation 
length. Therefore TCd = when this is not satisfied, which truncates the integrals over A and z in 
a manner that depends on the order of integration. When the above expression is satisfied, the ray 
in question is detectable, and = 1. 

To avoid internal reflection, the ray must meet the surface at an angle of incidence that exceeds 
the complement of the Cerenkov angle, where the Cerenkov angle is 6c = cos~^(l/nr) = 0.954 for 
Ur = 1.73. The transmission coefficient tu(P) goes to zero at that angle, which is when f3 = Tr/2 — 6c, 



and remains zero for all larger f3. Thus, the constraint in eq. (21 ) is already violated when internal 
reflection occurs, and there would be no formal need to include 'Hw,ct>' separately. However, grows 
so rapidly with angle that it is appreciably nonzero even just a few degrees from critical, and it 
quickly saturates near 0.7 for angles as small as 10 degrees from critical. That rapidly saturating 
behavior, along with the fact that appears in a logarithm so is only crucially important for 
weaker fields, motivates our choosing to model ty = 0.6 as a constant. This loses some accuracy 
in the result, but is much computationally simpler than following a function whose small-angle 
approximation breaks down completely for angles more than just 2 degrees from critical. Uniformly 
setting t|| = 0.6 implies that we no longer have zero transmission at the critical angle, so we now 
need to apply the total internal reflection constraint, Ti.w,(j,', separately. 

To determine if internal reflection is avoided, we need the cosine of the angle of incidence, 
which is given by the dot product of the surface normal h{'w, cp') and the ray direction 7(A, (f>, a). 
Here a is the angle the neutrino (and its shower) makes to the horizontal, where our convention 
is that a < is for downward neutrinos that are first encountering the Moon close to the shower 
of interest. Thus given the glancing incident neutrino angle a, which sets the orientation of the 
Cerenkov cone, and the electric field ray angles A and (p within that cone (where A = is at the 
Cerenkov angle and (p = is most directly toward the surface and least apt to internally reflect), 
we require 

h{w,<p')-j{A,^,a) > sindc , (22) 

which yields directly 

sin^c < sin(0c + A) cos 0sin fj cos 0' sin a — cos(0c + A) cos cr sin a 

— sin(0c + ^) sin(/)sin<Tsin^' + sin(6'c + A) cos </> cos o" cos a 

+ cos(^c + ^) sinacos 0' cos a . (23) 



Applied to the integrals in eq. (12), this constraint identifies a maximum azimuth (j) and truncates 
the (j) integral there, while the constraint in eq. (21) truncates the integral over depth z. 



3. Further approximations to reach a result in closed form 

The goal of this paper is to achieve closed-form expressions for the effective area A{E) as a 
function of neutrino energy, lunar regolith properties, and observing parameters. Since at this point 
we still have a cumbersome 6-dimensional integral to evaluate, we must avail ourselves of several 
additional approximations to obtain a tractable expression which yields the desired scaling laws. 



3.1. The near-surface emission approximation 



Since the integrand e""^" depends on z, the z integration is nontrivial, but this dependence is 
removed when the observable showers have to be so close to the surface that there is no appreciable 
change in the neutrino irradiation over the region where the radio signals are detectable. The 
requirement to make this approximation is 



« ^ 



(24) 



where and Ly are the electric-field attentuation length and neutrino mean-free-path respectively. 
This "near-surface emission" approximation allows the z-integral to be performed trivially. The 
expressions for Ly and are taken as (summing the cross sections for neutral and charged-current 
interactions, as they both generate similar hadronic showers, Gandhi et al. 1998 Reno|[2005 1 , 



0.07 



E 



"1/3 



P 



Rm VlO^OeFy \preg 

where preg = 1.8 g cm~'^ is the regolith mass density, and 



1 X 10" 



GHz 



-1 



(25) 



(26) 



(This result is known only to within about a factor of two ( [Gorham et al. 2004; Scholten et al. 
20061. Hence the approximation requires v ^ 5(£'/10^'^ eV)^/^ MHz for lunar rock with density 
p = 3 g cm^'^. This is satisfied over most of the range of neutrino energies and radio frequencies 
that concern us, except when the lowest frequencies < 100 MHz) are used to observe the highest 
energy (i?^10^^ eV) neutrinos. 



The surface-emission approximation allows us to replace Tu in e '^'^ by 



sin a 7i{a) 



(27) 
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for TC{a) again the Heaviside step function, so r^, = for q < (downward neutrinos) and = 
sino X 2Rm/Lv for a > (upward neutrinos). Here we insert the order-unity parameter ip = lA 
(appropriate for a neutrino spectrum ~ see Appendix) to account for neutrinos from higher 

energies being downgraded into the detection regime for energy E, which shghtly enhances the 
aperture. Note the absence of dependence on z, making the z integration a trivial matter of simply 
tracking its range of contribution. 



Carrying out the z integration subject to the selection rule in equation (21 ) and using 
s z 



sin(0c + A) cos ^ cos Q — cos(^c + A) sina ' 
the detection probability integral (eq. |12|) can be written 



(28) 



2 L /"oo 
P(E) = da cos a dA sin(6lc + A) 



oo 

roc /"OO 

X / # / #' / dwe-'^'Tizmax) e""" n^,^> m , (29) 

Jo Jo J-oo 

where Zmax is the maximum depth for which the emergent field exceeds the minimum detectable 
value, and we define 

T(/) = fHif) (30) 
for TC the Heaviside step function, to insure the integrand never contributes when it is negative. 



We by combining eqs. I\2lh and i28n 



IS? 

[sin(^c + A) cos (/> cos a — cos(0c + A) sin a] ( 1 — „ „ ) (31) 



Jo 



where fo is the dimensionless quantity 



In I |i ) . (32) 



Physically, fo is the ratio of the thickness of the Cerenkov cone corresponding to the electric field 
threshold to the 1/e full thickness (2Ao), as shown in Fig. [2| 



3.2. Small angle approximation 

We have now reduced the number of integrations to five, including two over surface roughness. 
One more integration, that over azimuthal angle around the Cerenkov cone (j), can be carried out 
surprisingly easily if we make the further approximation that the detectable radio rays hug tightly 
to the Cerenkov cone, such that A is small, and that the other angles we treat, (/>, a, and cr, are 
also small. The small-angle approximation is particularly valid for lower energy neutrinos near 



- 14 - 



1.2 




Emission Angle, 8^ = 54.7 

Fig. 2. — Cerenkov cone angular profile at E = 10^^ eV and v = 1 GHz. The dimensionless factor 
fo is the ratio of the cone thickness at threshold £min to the 1 /e width. 



the important GZK cutoff. It also applies to higher observing frequencies above 1 GHz, and may 
extend to lower frequencies as long as the highest energy neutrinos are not the focus. In this 
approximation, the maximum azimuthal angle that does not totally internally reflect, is given 



by equation (23) to be 



5i? 



2(A + waocoscj)' - a) ^^^^ 



y tan 9c 

which supplies the upper limit for the (j) integration. 

The small-angle approximation also permits us to use 

z^ax = sine, /2(^1 - , (34) 

and to replace sin by sin 6c = -y/l — and cos a by unity. This gives us all the expressions 
we need to simplify the evaluation of the integrals, recalling again that fo is influenced by t\\, 
and to within an expected accuracy of 10-20% we take t\\{(3) = 0.6, rather than its small-angle 
form, as the latter loses accuracy too rapidly for the angles we need to treat. We note that the 
primary influence of t|| is in determining the depth below the surface that will be visible, and 
since the angular contributions scale in a self-similar way as this depth is varied at different E, 
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the systematic errors introduced by our approach should appear primarily in the energy scale. 
Hence the transmission coefficient has its largest significance in the determination of the minimum 
detectable neutrino energy, and since this is already an important issue for GZK-type neutrinos, 
future work should attempt to clarify the reliability of the various and somewhat contradictory 
treatments of iii found in the literature. 



3.3. Evaluation of the azimuthal integral over the Cerenkov cone 

Approximating all remaining angles to lowest nonvanishing order simplifies the expressions 
dramatically, and allows us to carry out the (p and w integrations in closed form. Prior to those 
integrations, and subject to the approximations above, the detection probability becomes 

X dA T 1 - / da e-^^ / # mnd^R) , (35) 

J-oo \ Jo^oJ J -co Jo 

where in the surface-emission approximation we have 

r, = r{£] (30) 



and we have defined 



(37) 



as the angle through which upward (a > 0) neutrinos are successful at penetrating the lunar secant 
through rock with mass density p = 3 g cm^'^ up to the regolith layer (with ip = lA chosen in 
regard to a somewhat arbitrary neutrino spectrum with power-law index -2, see Appendix). 

To evaluate the (p integral, let us first let us define the angle e, which is the angle of incidence 
(from inside the Moon) to the normal to the surface, relative to the critical angle, so 

e = - Oc - p. (38) 

We replace /? by e because the latter is small and will be considered only to lowest nonvanishing 



order. To this order, we find eq. (11) becomes 



as) = "''^"'^J^'" ' (39) 



which demonstrates explicitly how the solid-angle magnification factor gets large as the critical angle 
is approached. This magnification effect is an important contributor to what would otherwise be a 
much smaller aperture. 
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The angle e is not one of the phase space integration variables, but is easily expressed in terms 
of those variables, in the small-angle limit: 



where in this limit we have 

2(A + wfJoCosc/)' - a) 

(^R = / ^ , • (41) 

ynf. — 1 



Now we encounter an interesting result, which follows from elementary application of eqs. (39) and 
(|40l): 

dcp i{e)'H{<PR) = -^n{<PR) (42) 

for any set of values for the other phase-space variables. In other words, in each equal-size phase- 
space bin over A, a, a, and (p' , the integral over (p yields either the same fixed numerical value, or it 
vanishes. Physically, this says that the solid- angle magnification effect perfectly compensates for the 
way internal reflection truncates the azimuthal window for ray escape. This can be no coincidence, 
and presumably would be more intuitively clear using some other choice of phase-space partition. 



3.4. Evaluation of the integral over incident neutrino angle 

The simple result from the (j) integration is especially helpful in carrying out the integral over 
a, because it adds no new dependence on a to the integrand. Thus the integrand remains simply 
g-TfCo^/o^o)^ which is amenable to closed-form integration. When a > (upward neutrinos), we have 
an exponential integral, and when a < (downward neutrinos), we have a trivial integrand. 



At this point it is convenient to introduce the scaled variables 

a 

V = — 
ao 

and 



(43) 



Jo'-^o 

and define the "roughness parameter" 

X = 7^ , (45) 

which characterizes how the window of acceptance of downward neutrinos can be expanded by 
surface tilts that avoid total internal reflection, and the "penetration parameter" 

Jo'-^o 

which characterizes the contribution of upward neutrinos that penetrate a long distance through 
the lunar rock to reach the shower point. Using these definitions, eq. (35) becomes 
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P{E) = Po{E)ac 



7r/2 



dcj)' 



dw e 



du (1 - n2) / dv e-^(^)W(i; - v^in) (47) 



where for brevity we define 



and 



Po{E) 



{nl - 1) (L 



L,. 



U WX COS (f)' 

y y 

The integration over v may be carried out explicitly, yielding 

rTT/2 



P{E) = Po{E) ao r 
Jo 



du {l-u'^)U+ + 



dw e 



W/Wo 



du {l-u^)U^ 



whereto = l/(x cos 0') and 



u WX cos ( 
U± = 1 + - ± 

y y 



+ e 



±(u/y — WX cos (f)' /wo) 



(48) 
(49) 



(50) 



(51) 



Although further progress can be made carrying out the elementary integrals over the Cerenkov 
cone angle u, the expressions become extremely unwieldy, and the two further integrals over the 
surface tilt parameters w and cj)' would require numerical integration anyway. So we pursue no 
further the formal integration of these expressions, and instead turn to calculating scaling laws 
in various asymptotic limits. We ultimately find that a convenient global approximation may be 
derived without significant loss of accuracy, by using the approximate method described next. 



3.5. Asymptotics and scaling laws 

Since the ultimate goal is to gain insight by deriving scaling laws and asymptotic limits, we 



seek an approximate evaluation of eq. (50) that is accurate in both limits of weak and strong 
roughness (x — > and x S> 1, respectively), which will then hopefully maintain approximate 
accuracy in between. We may accomplish this simply by evaluating the aperture in these two 
opposite extremes, keeping only the leading contributions in each regime, and simply adding them 
together for a global scaling law. Because the regimes of contribution are so physically separable, 
simply adding them produces results that agree with the full expressions within ~25%, and also 
yields physically insightful results. 



Taking x = produces closed- form results for the integral in eq. (50), but even those results 
are more complicated than necessary. Keeping only the dominant terms in the limit of either large 
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or small y, we find we achieve better than 25% accuracy with the simple asymptotic expression 

P{E)smooth = + y«<') • (52) 

Note the first term in the parentheses originates from downward neutrinos and the second from 
upward neutrinos, so we see the angular acceptance of downward neutrinos is controlled by the 
Cerenkov width when the surface is smooth, and the acceptance of upward neutrinos is due to the 
angle that can successfully penetrate the Moon. 

If we include roughness and consider the limit x 2> 1, we find that the upward neutrino detection 
rate is hardly affected, because upward detections are limited far moreso by penetration than by 
total internal reflection, giving them a completely different character from downward detections. 
However, the downward detection rate may be greatly enhanced by roughness, because downward 
neutrinos often produce Cerenkov cones that are almost completely internally reflected unless the 
surface encountered by the rays is favorably tilted, and we find in this limit 

. M_l)|l;3^,(_ii_.„) , (53) 

SO the role of /oAq for the downward neutrinos has been supplanted by ctq; rays avoid internal 
reflection not by pushing to the edge of the detectable Cerenkov cone, but rather by being lucky 
enough to encounter the surface where the tilt is highly favorable. 

Adding the smooth and rough limits thus yields the general expression, consistent with pre- 
existing inaccuracies owing to the idealizations and approximations already in place. Thus we find 
for the full aperture calculation 

A{E) - ^,i!|_ll^/3A„(M/,3 + + M/J , (54) 

where 

= foK (55) 
accounts for downward detections without help from roughness, 

accounts for downward detections assisted by roughness, and 

1 6 

*u = = 5.3 ao (57) 

accounts for the detection of upward neutrinos. 

The above results underline clearly the three basic angular scales of interest, the Cerenkov 
width Ac = /oAq, the surface roughness parameter ag, and the upward neutrino acceptance angle 
Oo, which with appropriate coefficients, map into the angular acceptance parameters ^ds^ ^dr^ 
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and The aperture for neutrino detection is simply dominated by whichever of these is largest, 
and note that the ao angle for upward neutrinos is benefited by an order-unity coefficient that 
is significantly larger than the others, owing to the much higher azimuthal angular acceptance 
when internal reflection is less of a problem. Nevertheless, the large scale of Uo (> 10° at > 1 
GHz) makes roughness an important contributor, especially at higher frequencies, whereas the /qAq 
parameter is largest at low frequencies. Upward neutrinos, on the other hand, prefer lower energy 
neutrinos that penetrate over a wider ao, but are often difficult to detect unless the instrument is 
extremely sensitive. 



4. Aperture Dependence on Lunar and Telescope Parameters 



The scaling laws provided by equation (54) may be viewed as the fundamental result of this 
paper. Each of the three terms has a simple physical interpretation that can be conceptualized in 
the following simple form: 

Pi = Po (58) 

for i = ds, dr, u, and where 

Po = ^^^^foAo (59) 

accounts for the inherent angular width of the Cerenkov emission as well as the fractional loss in de- 
tections owing to the overpenetration of neutrinos past the layers where radio rays may be detected. 
The terms represent effective acceptance angles for the "active" neutrinos, i.e., the neutrinos 
that can contribute to detections. The relative contribution of these terms to the total aperture 
is illustrated in Fig. [s] for a threshold field Smin = 0.01 /^V m^^MHz^^ and two characteristic 
frequencies, 150 MHz and 1.5 GHz. 



4.1. Dependence on Surface roughness 

It is apparent from Fig. [3] that the surface roughness contribution "^dr is the largest contributor 
to the total aperture at 1.5 GHz, whereas it is unimportant at 150 MHz. The roughness contribution 



can be very large at even higher frequencies e.g., the GLUE experiment at 2.2 GHz (Gorham 



et al. 20041, where the enhancement over a smooth lunar surface is a more than a factor of 3 at 
E = 10^^ eV. The contribution from surface roughness is important when it exceeds the sum of 
the contributions from the smooth surface and upward-going contributions For = 1 ZeV, this 
occurs at > 300 MHz with only a modest dependence on neutrino energy. The importance of 
surface roughness to the detection aperture at high frequencies has also been found in several Monte 



Carlo ray-tracing simulations of the escaping radiation from the lunar surface ( Gorham et al.pOOl 



Beresnyak 2003 James k, Protheroe 2009) 
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Fig. 3. — Total aperture versus neutrino energy at threshold field Smin = 0.01 /iV m~ MHz~ and 
two observing frequencies: 150 MHz (solid red line) and 1.5 GHz (solid blue line). The contribution 
from each of the three terms in equation ( |54[ ) is also shown. At both frequencies, downward-directed 
terms have the largest contribution except near the lower cutoff energy. For energies well above 
the cutoff energy, surface roughness has the largest contribution, except at low frequencies, where 
the smooth surface term already dominates without assistance from roughness. 



4.2. Dependence on minimum detectable electric field {£min) 

The minimum detectable electric field of a radio telescope with effective collecting area A^, 
system temperature Tgys, and bandwidth Az^, receiving linearly-polarized radiation, can be written 



(Gorham et al. 2004 1 



1 

Srmn = ( ^^^^'j ' V/m (60) 

where Nu is the minimum number of standard deviations needed to reject statistical noise pulses, kb 
is Boltzmann's constant, Zq = 377 0, is the impedance of free space, and Ur is the refractive index 
of the medium. This minimum electric field is plotted as a function effective telescope collecting 
area in Fig. |4j using nominal values of Nf^ = 4.0, Tsys = 120° K (assumed dominated by the Moon's 
contribution, limb pointing), and Ur = 1.73. 
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Fig. 4. — Minimum detectable electric field vs telescope effective area using eq. (60) and bandwidth 
Au =100 MHz (solid line) and 1 GHz (dashed line). See text for assumed values for other telescope 
parameters. 



Since the effective thickness of the Cerenkov cone Ac depends strongly on Smin, the total 
effective aperture will also depend strongly on Smin, as shown in Fig. |5j The left panel shows the 
aperture at an observing frequency 1.5 GHz, while the right panel is at 150 MHz. This plot shows 
the clear trade-off between collecting area and minimum detectable neutrino energy as the observed 
frequency changes: at high frequencies the minimum neutrino energy cutoff is lower, while at low 
frequencies the total aperture is significantly higher at a given fixed telescope sensitivity {Smin)- 

The aperture dependence as a function of both telescope sensitivity (Smin) and observing 
frequency is shown as a surface plot in Fig. |6j Panel (a) shows the total aperture for neutrino 
energies exceeding 10^^ eV, while panel (b) shows aperture values for E > 10^^ eV. The dark blue 
regions (zero aperture) correspond to the low-E cutoffs seen in Figs, [s] and [5} Note that for a fixed 
neutrino energy and telescope sensitivity {Smin), the maximum aperture results from choosing the 
lowest observing frequency that avoids the sharp cut-off region. 
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Neutrino Energy (eV) Neutrino Energy (eV) 

Fig. 5. — (a) Aperture vs. neutrino energy at observing frequency v = 1.5 GHz using telescopes 
with minimum detectable electric field fmm= 0.0001, 0.001, 0.01, and 0.1 /iV m-^MHz"^ . (6) 
Same as panel (a), but at observing frequency v = 150 MHz. Note the trade-off between collecting 
area and minimum detectable neutrino energy as frequency is changed. 



4.3. Aperture optimization 

One benefit of scaling laws is the ability to make optimization calculations using analytic 
derivatives. The aperture is energy sensitive, so to maximize the neutrino detection rate, we must 
make some assumption about the neutrino energy spectrum. We will consider two cases: a power- 
law distribution and a mono-energetic spectrum. 

First consider a canonical power-law neutrino spectrum that scales with neutrino energy as 
E~'^ . We wish to choose frequency v to optimize 

n oc j dE E-^P{E) . (61) 

Let us first assume that the surface roughness dominates the aperture, so we use ^dr, and note that 
this will be most appropriate at higher frequencies where the detectable Cerenkov cone is narrow. 
Then 

n ^ dE E-^^fl^o , (62) 



so 

3/2 



l + (i./2.3)i-23 ' + ^ 



(63) 



We can scale the frequency out of the energy integration using y = vE, and neglecting the weak 
dependence in the vj2.'i GHz term in the logarithm, yields that TZ scales nearly with Thus 
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Fig. 6. — (a) Aperture vs. minimum detectable electric field and observing frequency for neutrinos 
with energies E > 10^^ eV. (6) Same as panel (a), but for neutrino energies E > 10^^ eV. 



the optimal frequency to maximize aperture is as low as possible, at the expense of increasingly 
larger neutrino cut-off energy sensitivity and technical issues such as RF interference and ionospheric 
pulse dispersion. 



If, on the other hand, if we are targeting a particular energy Eq, we can focus TZ on just that 
Eo, and find 

E 



TZ (X u 



In 



+ k 



3/2 



(64) 



,l + (z^/2.3)i-23 

Evaluating the derivative of this expression, again neglecting the variation in the 1 + (z^/2.3) term, 
leads to the conclusion that the detection rate is optimized when the frequency is chosen such 
that fo = y^3/4 = 0.87. This means that for for detection of neutrinos at energy Eq, fo values 
substantially above unity are wasting aperture and would be better served by reducing v (and fo). 
Likewise, fo values substantially below unity are encountering the energy cutoff problem, so it is 
preferable to increase the observing frequency. 

At low frequencies surface roughness no longer dominates, and instead it is the broadly de- 
tectable Cerenkov cones that controls the aperture (^'^s) for higher energy neutrinos. Then we 
would find for an E~'^ power law 



7^ cx (IE E~^/^ 



In 



V E 
1 + (i//2.3)i-23 



-\ 2 



+ k 



(65) 
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which scales roughly like v and favors even more heavily the lower frequencies. However, if a 
particular energy Eq is the target, optimization occurs near /o = \/2/2> = 0.82. 

In summary, for power-law neutrino spectra with no high-energy cutoff, the detection rate 
is maximized at the lowest possible observing frequency, neglecting technical problems with RF 
interference and ionospheric dispersion. However, if a search is optimized for a particular neu- 
trino energy range, the aperture (and consequent event rate) is maximized when the dimensionless 
parameter /o 



0.8. From eq. (32) this condition can be written 



f 



1.9 



We can recast this condition to a more tractable form using eq. (18), 



■'opt 



1.23' 



43.2 GHz 



E 


-1 - 


ZeV 





f 

On 



(66) 



(67) 



_/xV m-i MHz" 

where Vopt is the observing frequency which maximizes aperture at neutrino energy E. 

Fig. [rj^a) shows the frequency which maximizes aperture for a given target neutrino energy for 
threshold sensitivities Emin = 0.001, 0.01, and 0.1 /iV m~^MHz~^ . Fig. [7|^b) plots the aperture 
as a function of frequency for a threshold sensitivity £min = 0.001 /uV m^^MHz"^ for three target 



energies near the GZK cutoff. As discussed in section 4.2, this sensitivity corresponds to much 
larger aperture (e.g., SKA) than current searches. Note that the optimal search frequency for GZK 
neutrinos depends strongly on the exact model for the cosmogenic neutrino energy spectrum, which 



varies with UHECR source model (e.g., Kalashev et al. 2002). 



4.4. Telescope requirements for GZK-neutrino searches 

In order to probe near GZK energies {E < lO^*' eV), the minimum detectable electric field 
must be £min ~ 10^^ /xV m^^MHz"^ , largely independent of frequency. From Fig. |4]we can see 
that even for very wide bandwidths, this requires an effective collecting area exceeding 50,000 m^, 
approximately the area of the Arecibo telescope. Since most detection schemes require coincidence 
on multiple telescopes to discriminate against accidental pulses, an array of multiple Arecibo-class 
telescopes are required. Current experiments using multiple 10-100 m diameter telescopes can only 
probe neutrino energies well above GZK, where more exotic neutrino production processes may 



be operating (e.g., Kalashev et al. 2002). We therefore conclude that only next-generation arrays, 
e.g., the full SKA (Ag ~ 10^ m^, £min ~ 10~^ /iV m~^MHz^^ ) will have sufficient collecting area 
to probe GZK energies. However, even the SKA will have an aperture of only a few km^ at the 
GZK cutoff Eqzk ~ 10^^'^ eV. Given an expected GZK neutrino flux Fj, ~1 km~^ yr~^ sr~^ (e.g., 
Gusev et al.||2006 ), the event rate might only be a few counts per observing month. Hence an even 



higher energy neutrino population, if it exists, might continue to prove easier to detect, even with 
SKA-class technology. 
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Neutrino Energy (eV) Frequency (GHz) 



Fig. 7. — (a) Optimal observing frequency vs. target neutrino energy for threshold £min = 0.001 
(dashed line), 0.01 (solid line ), and 0.1 (dotted line) m~^MHz~^ . (b) Aperture vs. observing 
frequency for Smin = 0.001 /uV m~^MHz~^ and neutrino energy E = 5- 10^9 eV (dashed line), lO^o 
eV (solid line), and 5 • 10^° eV (dotted line). 



4.5. Comparison with Monte Carlo simulations 



Finally, we address the question of how the present analytic calculations compare with previ- 
ously published Monte Carlo simulations (e.g., Gorham et al.|200l Beresnyak||2003 Williams|2004 



James & Protheroe 20091. In general, a direct comparison with each simulation is problematic. 



since the input physics model (neutrino and radio extinction lengths, Cerenkov cone width and 
peak electric field, lunar roughness, detailed regolith properties) is significantly different for each 
simulation. However, since we have largely adopted the input physics parameters from James &; 



Protheroe (20091, a direct comparison is warranted in this case. 



Table 1 shows a comparison of calculated apertures using the analytic approximation (equation 



54 ) compared with apertures reported by James & Protheroe (2009 1 in their Fig. 6a (no sub-regolith 
case, different pointings summed), at a fixed neutrino energy E = 1 Zev (10^^ eV) for convenience. 
For completeness, we also list the upper limit (90% confidence) to the commonly plotted quantity 
F{E) = E'^I(E), where I{E) is the differential neutrino flux, using the "model independent" 
expression ( Lehtinen et al.|[2004 ) 



F{E) < 



2.3 E 

tobsA{E) 



(68) 



where tobs is the total observing time. 
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Table 1 lists observational parameters, apertures, and flux limits for three searches: Parkes 



(Hankins et al. 1996 James et al. 2007a), GLUE (Gorham et al. 2004 1, and Kalyazin (Beresnyak 
et al. 2005). The table lists the relevant input parameters used in the calculation, which were 
derived from the literature: threshold electric field {Smin), center observing frequency {v), and 
fraction of the Moon's limb sampled (("). This fraction is difficult to estimate accurately since 
the searches often used pointings with roughly defined beam positions (e.g., 'limb', 'half-limb') 
and differing beamwidths for pulse coincidence schemes with multiple telescopes or frequencies. 
Nevertheless, we made our best estimate based on the experimental descriptions, and multiplied 
the calculated aperture by this fraction, since the calculation assumes 100% limb coverage. The 
columns in Table 1 are: (1) experiment name, (2) threshold electric field {Smin), (3) mid-observing 
frequency (z^), (4) total observing time (tobs), and (5) fractional lunar limb coverage {(). Columns 



6-7 are apertures calculated using eq. (54) (Agmj), and from Fig. 6a of James & Protheroe (2009) 
(Ajp). Columns 8-10 are upper limits to the quantity F{E) using aperture Agmj, Ajp, and directly 
from the original search papers. 



Inspection of Table 1 shows that the apertures calculated at £' = 1 ZeV by the analytic 



approximation agree very well with the Monte Carlo results of James Sz Protheroe (2009) for all 



three experiments. In Fig. [8] we show a comparison of the ratio of effective aperture to energy as a 



function of neutrino energy for the GLUE experiment using our analytic calculation (eqn. 54 ) with 



James & Protheroe (2009, regolith only). The agreement is quite good, although there appears 
to be a small systematic difference with energy scaling. It is difficult to assess whether or not 
any of the disagreement could be traced to weakness in the Monte Carlo simulations, as details of 
such calculations are not generally reported with sufficient completeness to reproduce the results 
in detail. Nevertheless, the overall conclusion is that the resulting estimates of neutrino fiux upper 



limits are in good agreement, and our results help confirm the assertion of James Sz Protheroe 



(2009) that the widely referenced GLUE upper limit (Gorham et al. 2004) should actually be an 



order of magnitude higher than previously reported. 

Table 1: Aperture, neutrino fiux limit comparison at neutrino energy E = 1 ZeV 



Experiment 


£min 


1/ 


tobs 


c 


Agmj 


Ajp 


lOg{Fgmj) 


log{Fjp) 


log ( For jg 




Vm"^ MHz-i 


GHz 


Hr 




(km^- 


-sr) 


(GeV 


— 2 —1 

cm s 


sr"-*^) 


Parkes (limb) 


0.013 


1.5 


2.0 


0.20 


2.1 


2.0 


-1.8 


-1.8 


-2.10 


GLUE (hmb) 


0.011 


2.2 


12C|^ 


0.17 


1.1 


1.0 


-3.4 


-3.2 


-4.2 


Kalyazin 


0.013 


2.3 


31 


0.11 


0.7 


0.6 


-2.5 


-2.3 


-2.2 



'James et al. 



^Upper limit: 



(2007b) - 



Hankins et al. 



Gorham et al. 



( 1996 ) does not give a fiux limit at 10 eV 



(2004) does not state how many hours were pointed on limb. 
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Fig. 



Effective aperture divided by neutrino energy versus energy for GLUE experiment 



(Gorham et al. 2004) calculated using equation 54 (dotted line) and Monte Carlo calculation of 



James Sz Protheroe (2009) for regolith only (open circles, from their Fig. 6a). 



5. Summary and Conclusions 

Accurate aperture calculations are complicated and laborious and historically have required 
Monte Carlo simulations. We have shown that an analytic calculation with a series of simplifying 
approximations can reproduce a result comparable to Monte Carlo simulations, and in addition 
generate simple scaling laws with a straightforward conceptual interpretation (eq. [54]). We find 
that it is crucial to account for surface roughness when apertures are intrinsically low, such as 
when using high frequencies from Earthlike distances, but when the aperture is intrinsically high, 
such as at low frequencies with high-energy neutrinos, or for lunar orbiters, then surface roughness 
is of lesser significance. We also find (Fig. [s]) that downward neutrinos significantly dominate 
over upward neutrinos for higher neutrino energies, and this is especially true for lower frequency 
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observations ^ 300 MHz). However, at the energies nearest to the GZK regime, which is 
perhaps of greatest cosmological importance, both upward and downward neutrinos may contribute 
comparably to the detection rate. 

These conclusions imply that the first detected GZK neutrino, if detected at high frequency 
to give a more favorable energy cutoff, may be a downward neutrino whose Cerenkov signal will 
have escaped by virtue of a lunar roughness feature, or it may with similar likelihood be an upward 
neutrino that will have made it through a long column of lunar rock and then created an upward 
Cerenkov flash that did not require assistance from lunar roughness. Alternatively, if the first GZK 
neutrino detection comes at lower frequency, it will likely be because a downward neutrino sent out 
a very weak radio signal that did not require lunar surface roughness to escape, but did require an 
extraordinarily sensitive radio instrument to detect. It is also possible that the first UHE neutrino 
will come from pion decay near the acceleration region of a UHE proton, and then it could come 
at an energy well beyond the GZK regime. In that case, the detection could come at low or high 
frequency, depending on whichever instrument first achieves the necessary sensitivity. 

Another important conclusion is that simple scaling laws, with all their extreme portability, 
may be used to optimize experimental design within a wide range of constraints. They are also useful 
for making instant comparisons between the apertures of past experiments using a standardized 
treatment of the relevant input parameters. In this way, the impact of the various inconsistencies 
between models, such as the detectable Cerenkov cone width, the treatment of transmission and 
solid-angle magnification, and the role of surface roughness, can be addressed without having to 
run Monte Carlo simulations with exactly the same parameters. Indeed, uncertainties in the details 
of such simulations make it difficult for us to resolve several inconsistencies between our results and 
others quoted in the literature, whether they be due to problems in the Monte Carlo simulations 
or in our own analytic approximations, especially our use of small- angle approximations and a 
constant transmission coefficient. For example, we have not been able to determine the source 
of the order of magnitude increase in the neutrino flux upper limit compared with the published 



GLUE upper limits (cf. Table 1). This discrepancy was also reported by James Sz Protheroe (20091 
using a Monte Carlo simulation. 

One particularly robust result we obtain is that to optimize the detection of neutrinos at a 
given energy, one should choose a frequency that will yield a value for the fo parameter of roughly 
0.8, which implies a ratio of the maximum to the minimum detectable field of about 2, for that 
energy. Inverting this implies that past experiments at given frequencies are best tailored to the 
neutrino energies for which the maximum field generated in the telescope by neutrinos of that 
energy is about twice the minimum detectable field for that instrument. 



We would like to acknowledge helpful discussions with Hallsie Reno, John Ralston, and Clancy 
James, who were especially helpful guides in areas of this effort outside our own expertise. 
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6. Appendix A: The effective neutrino extinction pathlength 



When UHE neutrinos strike the Moon, there is a roughly 1/3 chance (Reno 20051 that they 
will initiate via the neutral current interaction a hadronic cascade that will leave roughly 80% of 
the neutrino energy still in the highly forward-scattered neutrino. The remaining 2/3 of the time, 
a charged current interaction will destroy the neutrino, but still result in the same roughly 20% 
of the energy going into the hadronic cascades we wish to detect. Thus from the point of view of 
detecting the cascade, we do not care whether the interaction was via neutral or charged current. 
However, this is a relevant question when considering whether or not the neutrino survives the 
encounter and can initiate further cascades. 

Since the observed radio signals are from a narrow layer at the surface, it is highly unlikely 
that the same neutrino could be involved in multiple detectable showers. Nevertheless, some of 
the neutrinos we are hoping to detect (labeled the "upward" neutrinos) may have to make their 
way through a considerable amount of lunar material before arriving at their observable shower 
location. For these neutrinos, it does matter if they can survive prior neutral-current interactions, 
losing only ~ 20% of their energy each time. For example, a neutrino that suffers two neutral- 
current interactions and no charge-current interactions can arrive in the domain of interest with 
roughly 64% of its energy intact, so at the energy of interest, this just means that there can be a 
contribution from initially higher energy neutrinos. This "downgrading" effect is easily accounted 
for by increasing the effective neutrino extinction path over and above Ly{E), by a factor that 
depends on the steepness of the energy spectrum, thereby increasing the number of "upward" 
neutrinos that contribute to the aperture at any given energy. Once this has been accounted for, 
the actual reaction rates for a given neutrino population are still proportional to -LjJ^, so the effect 
is not the same as a global correction to Ly^ it is merely an effective reduction in neutrino extinction 
at energy E. 



This preprint was prepared witii tlie AAS M^rjX macros v5.2. 
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Here we derive in the simplest case the enhancement in the effective extinction pathlength, so 



we find the factor ip we use to multiply L^{E) in eqs. (27) and (37), for a power-law neutrino energy 
distribution, with flux per unit energy I{E) oc E~'^. The power-law nature of the problem imposes a 
scale invariance that allows us to seek a constant factor by which all effective extinction lengths are 
altered, where the extinction of I{E) must be modified by the appearance in I{E) of downgraded 
neutrinos from initially higher energies. This results in an effective extinction coefficient x{E), 
differing from the actual extinction coefficient such that I{s) = Ioe~^^. The net extinction 

over length ds must obey 

^ = -x{E)Ioe-^' + b- 1.077- 0.8 ■x{E)Ioe-^\ (69) 
as 



where b is the branching ratio for the survival of the neutrino, which here is b = 1/3 (Cooper- 
Sarkar k SaA^|2008[ ), 1.077 comes from assuming that xi^) oc E^^^ so the x for the downgraded 



neutrinos is slightly higher than x(-^)) ^-iid the 0.8 factor comes from the combination of the E~'^ 
neutrino power law (assumed only for simplicity here) and the fact that higher energy bins are 
squeezed into narrower energy bins when 20% of the neutrino energy is lost to the hadronic shower. 

It remains only to note that the above equation takes on the desired form 

when X = x/1-4. Hence, the extinction of neutrinos with I{E) oc E~'^ acts as though the extinction 
length was larger by a factor ^ = 1.4. Had we instead used a power law of -2.7 instead of -2, the 
factor would have been 1.3, but the factor could be higher for flatter spectra, such as in the rising 
hump of the GZK cutoff region. It is even possible for the neutrino flux at a given energy at the 
start of the GZK hump to experience a region of neutrino enhancement inside the Moon, but we 
do not deal with this possibility here because no suitably general assumptions about the shape 
of the neutrino spectrum can be applied at this time. If GZK neutrinos are detected, they will 
likely come at energies near the peak or the falling part of the GZK hump, so we are probably not 
underestimating the neutrino population that is being passed down from higher energies as they 
cross lunar rock. 



